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ABSTRACT 


“With growing Army-wide use of programmabl= calculators, 
a system is needed to support the programming and testing of 
calculator software. This thesis provides a Fortran IV pro- 
gram to enable an operations research analyst to nore effi- 
ciently write and document HP4ICV calculator programs. 
Optical kar code readable by «he HP41CV is generated by the 
program. Also given is an IBM EXEC II program which pro- 
vides an interactive programming environment including on- 
line, self contained instructions. To illustrate the use of 
the system and the quality of the finished bar code and cal- 
culator program listings, examples are given including sin- 
gle variable statistics and linear programming. A final 
example provides a set of short utility routines which 
illustrate how pregrams can be developed for use in a 
calculator read-only-memory. 
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I. ENTRODUCTION 


Por the Army to fight effectively in a resource scarce 
environment, the quantitative decision making techniques of 
operations research are important skills for Army staff 
officers. Staff officers are expected to be able to put 
nugbers in their estimates when briefing commanders. They 
are expected to be able to meaSur2= and evaluate complex 
operations and subordinate units. They are expected to be 
frugal managers of time and money. And above all, staff 
officers must be able to apply sound, quantified reasoning 
in planning how to win the air-land battle. 

The use of hand-held programmable calculators by Army 
staff officers has the potential for improving the use of 
quantitative decision making techniques throughout the Army. 
Paster and more accurate than paper and pencil, the calcula- 
tor is less expensive and more portable than larger comput- 
ers. Even when compared to the latest micre-computer 
systeas or to portable terminals used for distributed data 
processing, the hand-held programaable calculator offers 


advantages in cost, reliability, power consumption and emis- 


sion of electromagnetic radiation. Hand-held programmable 


ne in at a ee 


calculators have already been successfully used by soldiers 


in the field for applications in artillery fire direction, 
surveying, and navigation. In addition, large numbers of 
Aray officers own their own pocket calculators and routinely 
use them for staff planning and reporting functions. 

In January of 1981 the U. S. Army Command and General 
Staff College at Fort Leavenworth, Kansas selected a pro- 
gtammable calculator for the Combined Arms and Services 
Staff School (CAS 3.) Using both resident and non-resident 
instruction, this course is designed to teach all Army cap- 
tains staff techniques and procedures. As a significant t 
part of the curriculua, the students are introduced to sub- 
jects such as statistics and regression, decision theory, 
combat modeling and linear programming. Considering the 
large number of officers projected to attend this course in 
future years, this course represents the most widespread 
training in operations research techniques ever attempted by 
the Army. The decision to provide a sophisticated calcula-~ 
tor to these students on an experimental basis was made for 
two fundamental reasons. First, the availability of a cal- 
culator with immediate field utility should motivate the 


student to apply the quantitative techniques as compared to 


the student who would be forced to do all calculations by 


hand. Second, the power of the calculator permits classrooa 
discussion cf techniques such as linear programming ard 
regression which are very difficult and time consuming to 
perform manually. 

This thesis documents the author's work to suppor* the 
use of a calculator in the Combined Arms and Services Staff 
School. Initially, the intent was to produce a series of 
lesson materials incorporating the use of the calculator on 
a series of operations research topics which have immediate 
application for the Army iivision level staff officer. 
Instead, the work accomplished focused on the design and 
construction of a system to make the programming and testing 
of calculator programs easier and more efficient. Except 
for the introduction, this thesis is written for the person 
wishing to iaplement the prcegramming support systen 
described. The implementor must have a detailed knowledge 
of the instruction set and programming characteristics of 
the HPG1CVY calculator as described in Wickes [Ref. 1: pp. 
6-20]. For the eventual user of the system, as compared *o 
the implementor, the system itself provides on-line dccumen- 
tation on how to use the system and what commands and 
options are available. Figure 1 shows the command menu dis- 


played on the terminal screen by this interactive prograa; 


ee tencensne te 


Figure 2 gives a more detailed explanation of each of the 
commands; and Figure 3 displays the on line introductory 
material that is provided to new users of the system. Yor 
the user, a knowledge of the infermation contained in the 
calculator owner's handbook [Ref. 2] is sufficient to begin 
writing calculator programs using the support system 
described. 

The calculator selected by the Command and General Staff 
College, the Hewlett-Packard HP4ICV, typifies the state of 
the art in off-the-shelf calculator technology. While not 
without disadvantages, this calculator was selected because 
of its power and features which make it easier for Arny 
staff officers to use. Pirst and most important of these 
features is the ability of the calculator tc manipulate 
alphabetic characters in addition to numeric data. The cal- 
culator can display the name of a variable when input data 
is required or label output when the calculation is conm- 
pleted. With this feature, the calculator helps the user 
know what data to input or what action to take next. I¢ 
also helps alleviate the need for constant reference to 
printed instructions which are difficult to use under field 


conditions. 
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On-Line Introductory Material 


Pigure 3 


A second important feature. is the multiplicity of means 
by which programs can be entered into the calculator. Mag- 
hRetic cards, read only memory, and optical bar code are all 
available and each has advantages depending on the situa- 
tion. Fer the long term, read only memory offers the abil- 
ity to retain very large programs (in excess of 8000 bytes) 
and the simplest and most reliable means of entering pro- 
grams into the calculator under field conditions. For the 
short term, optical bar code offers the least expensive 
method of reproducing and distributing calculator scftware 
that has not been subject to extensive field «esting. In 
addition, as shown in this thesis, the optical bar code can 
provide an important link between a main-frame computer and 
the hand-held calculator. 

A third important feature of the HPYICV is its rela- 
tively large memory capacity as compared to programmable 
calculators such as the Texas Instruments TI-59. A large 
amount of memgory permits the solution of larger, often more 
realistic problems than could previously be solved ona 
hand-held device. A demonstration program given in this 
thesis for linear programming is an exaaple of an applica-- 
tion where the full memory capability of «he HPYICV is 


required to be able to solve realistic problems. 
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To take advantage of the calculator's unequalled ecc 
and portability, the operations research analyst is chal- 
lenged to overcome its limits of speed and memory capacity. 
The preparation of calculator software is as jifficuit, if 
not more so, than the preparation of software for larger 
computers. To accomplish the most possible with the hand- 
held device, the calculator programmer is often forced to 
write programs which are very difficult to comprehend when 
examined by other programmers. AS Dahl, Dijkstra and Hoare 
(Ref. 3: pp. 1-10] point out there are limits to human con- 
petence which interfere with the programming process. In 
the past, with less mature calculators which constrained the 
typical program to a few hundred program steps, these limits 
to human competence were neither as apparent nor as economi- 
cally important as they are with the HPY1CV. Accordingly, 

*t is not guvis toned that the average Army officer who uses 
the BPY4ICVY on real world problems which push the calculator 
to the limits of its capability would write their own pro- 
qrams. In particular, it was never intended chat the stu- 
dents in the Combined Aras and Services Staff School would 
be taught calculator programming. It is a tribute to the 


power of the device and the quality of the calculator soft- 


ware when a relatively inexperienced user can run coaplerx 


programs using little more than the digit entry keys and the 
run-stop key on the calculator. This does not mean that the 
user must not have a clear understanding of his problem or 


the solution technigue, but rather it means that the calcu- 


| 


lator should not require programming skill or extensive 
training prior to application. 

The growing complexity of calculator programs described 
above and the realization that calculator programs for Army 
field use are not programmed in the field, suggest the need 
for a system to support the development, distribution and 
maintenance of calculator software. An operations research 
analyst or other professional programmer must be able to 
more efficiently prepare calculator programs than by keying 
them into the hand~held device. By preparing the programs 
initially on a larger computer, such as the IBM 3033, the 
programmer can use the speed and storage capability of the 
larger machine to great advantage. In addition, the avail- 
ability of a full-screen video text editor speeds the pro- 
cess of program revision and maintenance. By providing a 
capability to integrate comments directly into the source 
code on the larger computer, program documentation is more 
easily provided. Essentially the idea is that a programmer 


would write the calculator program using a terminal 


a 


connected to a large computer. After the calculator prograa 
is entered into the large computer, a compiler progras 
running on the large computer would check the calculator 
program for errors and convert the mnemonic instructions 
into the “key codes" which are the numeric instructions 
actually executed by the calculator. Then an emulator 
program running on the large computer would take the numeric 
instructions from the compiler and execute the program--in 
effect making the large computer produce the same effects as 
the calculator only much faster and more efficiently for the 
peedeiabarcs Finally, when the program has been written and 
tested on the large computer, optical bar code is produced 
which allows for the economical distribution and use of the 
program in the field. [9 encourage the calculator progranm- 
mer to use the system described, this process should occur 
in an interactive programming environment in which the user 
can move from one step to another by issuing siaple commands 
such as those listed and described in Figure 2 and receive 


help or on line instruction whenever desired. Under «his 


proposed system, the advantages of both the larger computer 
and the hand-held calculator are used appropriately ina 
autually supporting manner. This thesis presents two of the 


components of this proposed system. First, an IBM EXEC II 


program is given which provides an interactive programming 


environment for users operating under IBM's Conversational 


Monitor System (CMS.) A short discussion of the design of 


this program and a complete copy of the sourc? code is 


contained in Appendix Dto this thesis. Secondly, a cross 


compiler written in IBM standard FORTRAN IV is provided for 


translating calculator mnemonic instructions into the key 


codes necessary for use by the emulator and also for the 


production of optical bar code. The term cross compiler 


rafers to the “act that the program runs on one machine (the 


larger compu*er) but compiles programs for another machine 
(the calculz* ir.) A discussion of the design of this 


program and a complete copy of the source code is contained 


in Appendix D. To make the program easier to understand and 


adapt to new requirements, it is aodularized into 24 


subroutines and is heavily coamented. 


To illustrate the use of “he system, two of the six 


example programs originally planned are provided in this 


thesis. Revised plans now call for the remaining four exaa- 


ple programs to be issued at a later date as Navil Postgrad- 


uate School technical reports. Because the reasons for the 


delay constitute some of the aost important lessons learned 


from this thesis research, Chapter 2 documents the process 
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with a «echnical discussion of the factors involved. The 
major conclusions described in Chapter 2 are sche need for a 
prioritized list of criteria with which to evaluate calcula- 
tor programs and the need for mor2 structure in the prograa- 


ming process. Chapter 2 is technically oriented and assumes 


meee e gm 


the reader is familiar with the concepts of structured 
prograaming. 
Each of the calculator program examples is described in 


a separate appendix in which the documentation listed in 


1. Progras Description 


2. Sample Probles 
3. User's Guide 
4. Source Code Listing with Comments 


5. Bar Code 


Figure 4: Components of Program Documentation 


Pigure 4 is provided. The first example on single variable 
statistics is documented in Appendix A and uses the calcula- 
tor in an area where calculators have long been used, but 


does so in a way that shows the unique capabilities of the 


’, 
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HPG1CV. A second example on linear programming is docu- 
mented in Appendix B and illustrates an area where calcula- 
tors have not received widespread application. Most 
calculator linear programs which have been published to date 
have been either incomplete algorithas or have been limited 
to very simple problems. 

A third example, which by its nature does not conform to 
the documentation standards outlined above, describes a set 
of utility routines which could be distributed in read only 
memory. Programs for read only memory have different char- 
acteristics from other calculator programs and Appendix c is 


provided to illustrate some of these differences. 
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A. CHAPTER OVERVIEW 

This Chapter examines calculator programming within the 
context of the author's experience in preparing HP4ICV pro~ 
grams in support of the Combined Arms and Services Staff 
School. With the advanced capabilities and features of the 
HPWICY, it was hoped that a complete package of software 
could be prepared quickly. To document why this did not 
occur, this chapter will examine strengths and weaknesses of 
the calculator in relationship to a collection of techniques 
raferred to in computer science as structured programming. 
For the reader unfamiliar with this term, the previously 
cited work by Dahl, Dijkstra, and Hoare [Ref. 3] is 
recommended. This chapter is technically oriented and does 
assume familiaritiy with structured programming concepts. 

When prograarging calculator programs for personal use, 
most programmers, including the author, do noc find the «ask 
difficult. Programa@ing 2 hand-held calculator with the 
capabilities and features of the HP4ICV can be a rewarding 
experience. It is rewarding to master the algorithm of an 


operations research technique on a hand-held device. The 


aducational value in programming the calculator has been 
recognized by many educators, including damaiag [Ref. 4: 

pp. 2-3] and Weir (Ref. 5: pp. xii-xiiij. Providing a pro- 
gram for general distribution which makes optinum use of the 
calculator is quite a different situation. I+ was the 
author's axperience that programs, which gave correct 
answers when used by the author, often had to be compietely 
re-written several times before being acceptable. This 
problem became more acut2? as the size of the programs grew 
beyond 400 program steps, for at that size it became 
increasingly difficult to modify programs without affecting 
the total design. The major conclusions described in this 
chapter are the need for a prioritized list of criteria with 
which to evaluate calculator programs and the need for more 


structure in the programming process. 


B. STRUCTURED PROGRAMMING WITH THE HP4ItCYV 
1. The Need for structure 
To increase the efficiency of the programming pro-~ 
cess, a collection of techniques known as structured pro- 
gramming has received widespread attention in the computer 
science comaunity. While there is no one definition of 
structured prograaming, it does require three essential 


characteristics. First, there must be a logical structure 
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to the program which reflects the nature of the problem to 


be solved and any constraints imposed upon the soluticn. 
Second, the systematic process of stepwise refinamert is 
used to limit the complexity of program segments. Third, 
the programming language must reflect «he logical structure 
of the program and assist in stepwise refinement. These 
three characteristics represent not so much a detailed rec- 
ipe for program development as they do a philosophy of how 
programs can be more efficiently written. It was with this 
philosophy in mind, that a calculator programming support 
system was proposed which could take into account the 
strengths and weaknesses of the calculator; balance the 
structured prograsming philosophy with the other criteria 
listed below; and thereby solve the problems encountered in 
writing calculator software for the Compined Arms and 
Services Staff School. 
2. Fupdamental Limitations of Calculators 

Writing programs to solve complex problems on a 
hand-held calculator is difficult both because of inherent 
limitations in the calculation speed and memory capacity of 
the machine and also she inability of the calculator's 
native programming language to directly support structured 


progragaming constructs. In many respects, the task is 
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Similar to writing assembly level language programs fer 
larger computers. Calculator progranuming features a 
powerful instruction set including advanced mathematical 
functions but lacks any ability to refer +o variables by 
name instead of storage address. Like assembly language, 
the calculator's programming language consists of short 
mnemonic instructions typically followed by the storage 
location of the data to which the operation is to be 
applied. While a large amount of computer programming is 
still done in assembly language, it is generally accepted 
that programming in a higher level language such as PORTRAN 
is preferable. Programs written in an assembly language 
take more time to write and are not as easily changed as 
higher level language programs. Also, because? they depend 
on the instruction set of a particular machine, they can not 
be easily transfered from one computer to another. These 
same disadvantages apply to calculator programming. In 
addition, because the hand-held device does m>% have tne 
speed and memory capability of the larger machina, the 
calculator programmer must be even more mindful of the need 


to optimize his program to save program steps and execution 


tiae. 


3. Nodular Design 


The HP4ICV supports structured programming as well 
or better than any other hand-held calculator. As described 
in the owner's manual (Ref. 2: pp. 177-196], the machine 
primitive instruction XEQ encourages the conscruction of 
modular programs using calculator subroutines. Each subrou- 
tine can be a self-contained unit capable of being written 
and tested independently and used by multiple programs. 

This modularity is most strongly 2ncouraged wher routines in 
read only memory are used, for then the application progran- 
mer can significantly reduce the number of program steps in 
his own program. This modularity, however, is rot complete, 
since all variables are globally referenced and can be 
changed daliberately or inadvertently by any subroutine. 
This problem is no more apparent than with the use of read 
only memory, since one of the most limiting factors in using 
the read only memory programs as subroutines is conflict in 
the use of common registers. Also, unlike the modularity 
required in truly structured programs, there is no 
restriction limiting a subroutine to a single entry and a 
single exit point. In structured programs, such limits on 
entry and exit serve to define the fundamental building 


blocks by which stepwise refinement is made possible. With 
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the calculator, however, multiple entry and exit points are 
most useful for allowing a common routine to handle a 
duplicity of problem conditions. In this thesis, for 
example, programs are given for which two standard entry 
points are provided. One entry point uses an alpha-numeric 
label and an audio prompt to speed data entry, while a 
second entry point uses the alpha-numeric label but 
Suppresses the audio tone. After data entry, the value 
entered is displayed, and the user is required to verify the 
accuracy of the data entered. By using the same subroutine 
with different entry points, memory space is saved overall 
at the sacrifice of the structured programming philosophy. 

4. Control of Program Flow 

A basic deficiency prohibiting che HP4YICV from 

directly supporting structured programming is the way in 
which program flow is controlled. Programming languages 
which support structured programming typically have instruc- 
tion constructs such as WHILE-~ENDWHILE, REPEAT--UNTIL, or 
LOOP--QUIT-~ENDLOOP which make programming loops clea# and 
concise. Constructs such as IP-~THEN--ELSEIF--ELSZ-- ENDIF 
and tr> CASE statement make the evaluation of conditional 


expressions efficient and relatively error free. Also, 


structured programming languages typically discourage «he 
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use of GOTO unconditional transfers because they lead to 
confusing code. In contrast, the HPYICV programmer must 
woite his own looping constructs and his own corditioral 
evaluation constructs using machine primitive instructions 
which somewhat obscure the program's basic objective and 
flow of control. In addition, it is difficult to avoid dis- 
turbing pending operations in the stack registers when a 
conditional statement must be evaluated. AS can be seen by 
the short programa shown in Figure 1, the notation of the 
programming language does not permit structured program 
flow. 
5S. Clarification of Program Structure 

Because no calculator, including the HP41CV, sup- 
ports named variables, the use of comments as an integral 
part of the calculator program is vital if the logical 
structure of the program is to be made clear as required by 
structured programming. Comments should provide the vari- 
able names when storing and recalling data; they should pro- 
vide clarification of program flow; and they should mark 
subroutine boundaries and entry and exit points to make it 
easier to identify segments of the programa. With the 
HP@ICV's stack oriented architecture, it is also frequently 


useful to display the names of the contents of each of the 


Given the number n in the x-register, this program 
fragaent will sum the data values stored in memory 


locations 1 through n. 
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iastzuction Comsent 

LBL "SUM To execute press "XEQ SUM". 
1E3 

LG 

1 

| + Establishes a loop counter. 
0 Clears x and pushes loop 

| LBL 00 counter into y. 
RCL IND Y Recall the next data value. 
+ Accumulate the sun. 
Isc Y Increment the loop counter. 
6TO 00 If more data remains, branch; 
RIN else, quit and display sua. 
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Figure 5: Example Program to Add n Numbers 
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In Appendix C on common subroutines with 


stack registers. 


read only memory application, a shell sort [R2f. 6: opp. 
84-95] routine is given which employs che technique of using 
comments to display the names of the variables on the stack 
register. 
6. Data Types and Indirect Addressing 

Calculator programs represent more than a sequence 
of keystrokes; they also represent the manipulation and 
transformation of data. For maxinum efficiency, the manipu- 
lation of data should be structured so as to prevent common 
programming errors. For this reason, most computer lan- 7 
guages which directly support structured programming enforce ~ 


data type correspondence between data and operations. Pre- 


quently the formal declaration and initialization of vari- 


The HPYICV handles two types of 


ables is also required. 


data--real numbers and alphanumeric characters. While no 


formal declaration of variables is required, «ype checking 


is done automatically and is transparent to the user. Any 


attempt to perform an arithmetic operation on alpha-numeric 


data will result in the message "ALPHA DATA" and the program 


will halt. 


Because there is no formal declaration of variables, 


the programmer writing programas for the HP4SICV gust use 


extreme caution in managing his data set and insuring that 
the numbers stored and recalled by the calculatcr program 
are in fact the data elements desired. A typical example of 
an improper data reference occurs when a program is using 
indirect addressing and attempts to store or recall data 
frem a non-existent data register. This programming error 
is so common that a special error message "NONEXISTENT" is 
provided by the calculator when this error is detected. 
Indirect addressing is an important feature which gives the 
calculator a considerable amount of power and flexibility, 
but also represents an additional responsibility for the 
programmer to explicitly control. On the HPWICV all “indi- 
rect addressing calculations must be specifically provided 
by the application program--there are no vector or array 
data types such as usually found with higher level lan- 
guages. In an attempt to make indirect addressing more 
transparent to the programmer, an experimental subroutine 
was prepared to recall an arbitrary element of a matrix 
stored as a two dimensional array. This subroutine, which 
is shown in Pigure 2, was used in a simultaneous differen- 
tial equation combat aodel and the results evaluated. It i 
accomplished the task, but slowed the execution of the pro- 


gtam considerably (resulting in an overhead of 10.5 seconds 
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of axtra execution time for every 100 subroutine calls) and 
did not significantly improve the size or legibilixcy cf the 
application program. Accordingly, this technique is not 

recommended and indirect addressing remains a task that must 


be treated explicitly by the application programmer. 


C. ADDITIONAL CRITERIA FOR PROGRAM EVALUATION 

Calculator programming in many respects resembles a 
multi-criteria decision problem. On the surface the 
criteria for program effectiveness are quite straight 
forward--the program must yield the correct answer, run 
quickly, require the fewest possible memory registers and be 
user friendly. Unfortunately, these objectives often 
conflict and can not always be simultaneously achieved. In 
particular, the principles of structured programming are 
often in conflict with the desire to reduce the size of 
programs and increase their execution speed. It is also 
true that the objectives of structured programming concern 
the process of writing programs, whereas the additional 
criteria listed concern the final program product itself and 
ara therefore logically considered separately. Attempting 
to achieve all criteria at once can lead to failure, and 
some tradeoffs gust be considered to evaluate programs and 


guide program development. The following criteria represent 
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Entry to this routine assumes the x register contains | 
the column number and the y register contains the row 
number. The base address must be stored in ROW and 


the dimension of the matrix must be stored in ROS. 
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1 LBL "RCUM 
2 RCL O04 (BASE ADDRESS REGISTER | 
3 + (ADD BASE TO COLUMN NUMBER | 
4 X<>Y (RECALL THE ROW NUMBER 

5 1 | 
bs | 
7 RCL OS (DIMENSION OF THE MATRIX | 
8 * { 
9 + (ADDRESS IS NOW IN K REG | 
10 RCL IND X (RECALL THE DATA DESIRED 

11 RIN 

12 END 


Pigure 6: Program to Recall an Element of a Matrix 


34 


eee Lo 


“Lessons learned" in develorfing application programs as 
examples for this thesis. 
1. User Exriendliness 

User friendly programs consider the application 
environment and do not task th2 user to be all knowing or 
without error in entering data. While individuals differ 
greatly with experience, the average user will make frequent 
errors in entering data with the hand-held calculator's 
small keyboard. In talking with officers who had used the 
TI-59 calculator in the field for fire direction, it was 
discovered that most preferred to use the printer with the 
calculator kecause it allowed data to be checked after 
entry. This was in spite of the fact that the printer and 
calculator combination is more costly, less portable and 
less suitable for use in the field than the calculator 
alone. In short, user friendliness was more important than 
these other criteria. For this reason, it should be manda-~ 
tory that any calculator prograas intended for Army use in 
the field must allow the verification of data after entry. 
Because the use of the printer obviates many of the advan- 
tages of the hand-held calculator, the printer should not be 
reguired for this verification. Ine of the considerable 


advantages of the HP4ICVY is that the large amount of prograa 
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memory makes it possible to store the input values and 


perform this verification. However, programs written with 
this critaria in mind may not appear be most 2fficient to 
the casual observer. 

Another important aspect of user friendliness is 
liaiting the complexity of the calculator and the actions 
required «=o get results. The typical Army officer has lit- 
tle appreciation for the multitude of scientific and mathe- 
matical functions labeling the keys of the HPYICV. Yet the 
common programming practice of using the top two rows of 
calculator keys to indicate the identity of a variabie 
either upon input or output increases confusion over the use 
of the function keys. This works as follows: When local 
alphabetic labels are used in a program to represent entr 
points by which a user indicates the identity of an input 
variable or requests a particular output variable, then «he 
first two rows of keys on the HPU1ICV become subroutine exe- 
cution keys pointing to these local labels when the calcula- 
tor is in user mode. This feature was very important on the 
HP67 and TI-59 where the lack of alpha-numeric capability 
required «his method of program execution in order to most 


easily determine the identity of the input or output value, 


but it is less important on the HP4YICY. It is almost always 
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true that a program which requires the use of local labels 
is harder to use, and requires more frequent reference to 
the user instructions than a program which uses only «he 
run-stop key and properly prompts the user and labels output 
values. 
2. Execution Speed 

The second most important criteria for a calculator 
program is that it must yield results relatively quickly. 
In preparing example programs for this thesis, this point 
became very clear when testing two particular pregrams. One | 
program, a Simultaneous differential equation combat model, 
required in excess of 150 data values in order to yield 
results. It should be noted that it was only with the 
introduction of the HP4ICV that it became feasible to con- 
sider such large problems on a hand-held device. To accomo- 
date the size of the model, the program was written so as to 
economize on program steps at the expense of increased exe- 
cution time. It became immediately obvious upon initial 
testing that this had been the wrong priority--for users of 
the program were not impressed with either the use of the 
calculator or the utility of the combat model. If such user 
acceptance is not present, then the calculator program will 


remain unused, no matter how elegant tha design to conserve 
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memory. In contrast, the linear programming =xample given 
in Appendix B was written so aS to emphasize speed even i 
it meant including code redundancy. This program has been 
well received in part because it is so much faster than 
paper and pencil methods. 

The easiest and most 2ffective technique that is 
useful in increasing speed is to decrease the number cf pro- 
gram steps that the calculator must process inside program 
locps. For example, if two different program options 
require similar but slightly different actions within a pro- 
gram loop, it is tempting to insert a program flag check and 
branching instructions within a loop so as to use the same 
loop for both conditions. But this means that the calcula- 
tor must test the flag and branch inside the loop evan 
though the program is probably shorter overall.! 

Instead, if the application permits, the memory capacity 
of the HPYTCY can be used to best advantage by testing the 
flag once and then providing separate program loops for the 
two conditions. Again, this does not appear 2legant +o the 
casual observer, but it may result in a more successful 


program overall. This principle was discovered while 


1 Bran bing is reguired when ee fl tests either set or 
Clear if madre than one instruction i ,creguired to account 
for the differences in the two conditions. 
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programming the single variable statistics program given in 
Appendix aA. Initially, this program used a common loop for 
all data input and output operations, including reviewing 
the input data and making individual corrections. By 
providing a separate, somewhat redundant loop for data 
correction, the time required to input data points was 


reduced. 


D. A PROGRAMMING SUPPORT SYSTEM 

Considering the structured programming philosophy dis- 
cussed above in paragraph B and the additional criteria for 
evaluating programs listed in paragraph C, it becomes imme- 
diately obvious that programming with the calculator alone 
will never meet even a majority of these objectives. It 
must be recognized that the problem under consideration is 
not how the average person who owns a calculator should pro- 
ceed to program it for his own personal use, but rather how 
the Army can best provide the most cost-effective computa- 
tional resource for field use. For these reasons, a 
comprehensive programming support system is required. The 
programping support system outlined here will consider only 
the requirement for cost-effective preparation and 
maintenance of the calculator programs and not the breader 


issues of distribution and logistic support for the entire 
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calculator system to include hardware, training materials 
and printed references. 

1. A Cross Compiler and Bar Code Generator 

The first requirement for an operational support 

package is to free the programmer from the limitations of 
the hand-held calculator itself. Even with the printer and 
other peripherals, the calculator is no match for the larger 
machine when large programs must be examined or edited. In 
addition, the calculator is not currently capable of aceauce 
ing its own optical bar code as required for aconomic repro- 
duction and distribution of the software. Accordingly, a 
cross compiler for the HP41CV was listed as the first 
reguirement of the programming support system. Such a cross 
compiler has been written and is the major outcome of this 
thesis effort. This cross compiler accepts an HP41CV pro- 
gram written in the language of the calculator and returns 
the finished bar code as output. Any valid HP4Y1CVY program 
will be processed without need for modification by the cross 
compiler. In addition to the basic language of the calcula- 
tor, the user is allowed to inject comments directly into 
the source code with the use of the left parenthesis as a 
comment indicator mark. The ability to make comments 


directly in the source code sakes the calculator programs 
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more legible and more easily modified at a later date or by 
another programmer. Often, well placed comments can make up 
for a lack of structure in the program itself as far as 
legibility and maintainability are concerned. Having the 
comments directly in the source code facilitates their use 
and helps insure that they are as up to date as the progran. 
For the average programmer, use of unmodified HP4Y1ICVY source 
code augmented with a comment indicator will represent the 
most common use of the cross compiler. The cross compiler 
is described in more detail in Appendix D including a 
complate listing of the source code. { 
2. 4 Calculator Emulator Z 

After the calculator source code has been processed 
by the cross compiler, a need exists to be able to run «he 
program without the wait for the generation of bar code. In 
addition, for the future development of read only memories 
for the calculator, an emulator program is required because 
the calculator itself can store only up to 2000 instructions 
in active random access memory. The read only memory can 
store up to four times this amount. Thus, the calculator by 
itself may not be capable of testing extremely large pro- 
gvams or Renner with large amounts of constant data also 


stored in the read only aemory. Although an emulator was 
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not written for this thesis, the design of the cross 
compiler reflects the need for such a program. For example, 
the cross compiler generates an intermediate array of deci- 
mal integers which repesent the machine language of the 
HP&1CV prior to conversion to binary. It was intended that 
these decimal integers could be used without modification or 
further translation within a FORTRAN computed goto state- 
ment. Thus, with the difficult translation, instruction 


parsing and syntax recognition already performed by the 
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cross compiler routines, the emulator could consist of one 
large FORTRAN loop wherein a decimal integer was addressed t 


in the instruction array by a program pointer variable. The 
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integer is then immediately sent to a computed goto state- 
ment which would branch to the appropriate line of FORTRAN 
code which would simulate the referenced instruction, 
including updating the stack and the program pointer as 
appropriate. 
3. A Higher Level Language Compiler 
The final component in the calculator programming 


support system would be a program that would translate a 
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higher level language such as PASCAL into HP4ICV language 
which could then be sent to the cross compiler for verifica- 


tion and generation of the bar code and intermediate af 
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calculator language listings. It is the higher level 
language compiler that would most directly make up for ‘the 
weakness of the calculatsc in supporting structured prograa- 
ming. It would be able to increase the modularity of pro- 
grams, provide for named variables, make indirect addressing 
transparent and provide structured statements such as 
WHILE-~ENDWHILE and IP--THEN--ELSE. Again, the design of 
the cross compiler anticipates this requirement and provides 
a consideratle number of subroutines that would also be 
required by a higher level language compiler. These subrou- 
tines include a complete set of string functions for manipu- 
lating character data in FORTRAN and arn instruction parser. 
Because it was envisioned that the higher level language 
compiler would also be able to process statements entered 
directly as HP41CV instructions, the cross compiler is con- 
structed so that the routine which compiles individual lines 
of HP41CV source code could be called as a subroutine by che 
higher level language compiler. Thus, all three major coa- 
ponents of the proposed calculator programming support 


system would work together efficiently. 
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APPENDIX A 
SINGLE VARIABLE STATISTICS EXASPLE 


INTRODUCTION: 

Calculating single variable statistics is one of the 
most frequently used applications of programmable calcula- 
tors. Army division level staff officers use single vari- 
able statistics to summarize and describe data for command 
briefings and periodic reports. The text by Mendenhall, 
Scheaffer and Wackerly {Ref. 73 pp. 3-13] is recommended as 
an introduction to the statistical measures calculated by 
the program given in this appendix. This program 


automatically calculates: 

* Mean and Median 

e Sample Standard Deviation 

e Sua of the Squared Deviations about the Mean 
e coefficients of Skewness and Kurtosis 

e MSiniaua, Saximua and Range 


e Histogram Cell Frequencies 
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A single variable statistics program has been given as 
an exaaple because of its immediate utility to the staff 
officer and to illustrate several features of the HP4ICV 
which make it a superior device for Army field use. The 
most important of the these features is alphanumeric prompt- 
ing for input data values. The program given in this 
appendix provides an alphanumeric prompt for every input and 
output value and requires only the digit entry keys and 
run/stop key for data entry. Another important feature of 
the HP4ICV used by this program is its large memory 
capacity. This program retains up to 219 data points in the 
calculator's memory to allow the user to review the input 
data and make corrections during data entry. The large 
amount of memory allows the calculator to sort the data and 
calculate the order statistics including the mininua, 
maximum and median. Calculation of the median is a feature 
of this program which distinguishes it from other calculator 
statistics programs. In addition, without having to 
re-enter the data, the histogram may be calculated with a 


varying number of cells or a varying cell width. 


PROGRAM DESCRIPTION: 
The single variable statistics program has entry points 


for two different techniques of data input. The fastest 
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method, which provides both an alphanumeric prompt and an 


audio tone to speed data entry, may be called by execution 


| of the program from entry point "STAT1." A slower method, 
| which provides greater accuracy and suppresses the audio 
| tone for classroom use, may be called by execution of the 


program from entry point "$1." When called from "S1," the 


program requires the verification of each data point after 
entry. The sequence of actions is as follows: 


1. The calculator displays an alphanumeric prompt. As 
an example, "X1?" is the prompt for the first point. 


2. The user enters the data value with the digit entry 
keys and presses the run/stop key. 


3. The calculator displays the data entered with a label f 
derived from the alphanumeric prompt. For example 
"123.1415" is a typical calculator response. this 
display is prompting the user to verify the correct- 
ness of the data displayed. 


4. If the value is correct, then the user simply presses 
the run/stop key and the calculator advances to the 
next point. 

5. Ifthe value is erroneous, the user enters the correct 
value with the digit entry keys and then presses the 
Tun/stop key. Thén the calculator will again repeat 
step 3 and ask the user to verify the data value. 

This process will continue until the user makes no 
modification to the data value. 


To trun the program from 2ither entry point the user may use 
the XEQ key, or assign the entry point label ("STAT1" or 
"S1%) to a key and execute it by pressing that key in «he 
USER mode. Further instructions on running programs and 
making key assignments are contained in the calculator 


owner's manual (Ref. 2: pp 114-116). 
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In addition to the two initial entry points described 
above, several other alphabetic labels provid? the user with 
functions that are called outside the normal sequence of 
program execution. Label "SR" provides the user with the 
capability to review the data stored in calculator memory, 
aither before or after the data has been sorted. When used 
before the sort, the "SR" function is most useful in verify- 
ing the entire data set at one time. If used for this pur- 
pose, it should be called after all of the data has been 
entered and the mean of the data set is displayed with the 
"XBAR" label. If flag 21, the printer enable flag, is set 
"on" during this data review, then the calculator will stop 
as each point is displayed and the user may make corrections 
in the same manner as described above for the point-by-point 
verification associated with the "S1" entry point. When 
used after the sort, the "SR" function is most useful for 
displaying the order statistics for the data set. If used 
for this purpose, it should be called after the histograrg is 
output--when the "CMD" prompt is displayed. If the user 
presses run/stop after the "CMD" prompt, the order 


statistics will automatically be displayed. 


The design of the prograsa, especially the data entry 


loop, reflects the need for calculation speed. Code 


redundancy exists at several points in order to reduce “he 
need for extra flags, labels and goto statements which would 
slow execution during data entry. In spite of this need for 
speed, the summary totals needed for calculation of mean, 
standard deviation, skewness and kurtosis are accumulated 
during data entry. This is done so that these summary 
statistics are available with little or no wait following 
data entry. 

A complete listing of the program registers and flags 
used by this program is shown at the end of the program 


listing. 


SAMPLE PROBLEM: 
In order to establish a training standard for an obsta- 
cle coursa, a division assistant 33 randomly selects 10 sol- 


diers and records the time it takes each to complete the 


course. The following times in minutes were recorded: 


Determine the summary statistics and cell frequencies 


necessary to plot a histograa of this data. 


SOLOTION: 


1, First, set the size of the calculator's data nemory 
large enough to retain the data values, This requires at 
least 16 registers plus 1 for each data point, or a total of 
26 in this example. Alternatively, the size of data memory 
may be set arbitrarily large, up to a maximum of 235 pro- 
vided the user has no other programs in the calculator he 


wishes to retain. Por this example press: 
XEQ ALPHA SIZE ALPHA 26 


2. To call the program, determine the appropriate method of 
data entry and select the corresponding entry point. 
Press: 


XEQ ALPHA STATI ALPHA (quick entry) 
XEQ ALPHA $1 ALPHA (classroom use) 


3. The calculator will respond with the prompt "N?" asking 


for the number of data points. Press: 
10 R/S 


4. The calculator will respond with the prompt "X1?" asking 


-” 


for the first data point. Press: 


2.1 R/S 


Steele cea eel td 


NE eee kot om 


5. If you called the program via "$1" the calculator will 
respond with "X1=2.100" asking for verification that the 
first point is correct. If not correct enter the correct 


value, else press run/stop. 


6. The calculator will continue in the same way as steps 4 
and 5 for the remaining data points until all the data has 
been entered. If at any time you discover that you have 


made an error in data entry for any point, pr2ss: 


XEQ ALPHA sc ALPHA 


The calculator will respond with the prompt "POINT?" asking 
for the number of the point in error. For example, if point 


number 5 were in error, you would then press: 


5 R/S 


Assuming you had just input a5 as the point in error, the 
calculator would then respond with the prompt "X52?" asking 
for the correct value of point 5. Respond with the correct 
value and press run/stop. The calculator will then go back 
to the place in the data entry sequence where it left off or 


it will go to the calculation of the summary statistics if 


data entry was previously completed. 


7. when data entry has been completed, *he calculator will 


respond with the mean of the data sample labeled as follows: 
XBAR=2. 470 | 


At this point, you have the option of reviewing the entire 
data set or continuing to calculate the remainder of the 


statistics. To review the entire data set, press: 
XEQ ALPHA SR ALPHA 


Note that if flag 21 is set on (press SPF 21), the calcula- 
tor will stop after each data point is displayed, permit+ing 
you to change any value simply by entering the new value and ! 


pressing run/stop. 


8. After the mean is displayed with the "XBAR" label, if 
you simply press the run/stop key, the calculator will cal- 
culate the following statistics with the label shown: After 


each press R/S. 


Display Meaning 
$8 QD20.521 Sum of Squared Deviations 
About the Mean 
ae 0.24 Sample Standard Deviation 
EW= 33°} 70 Skewness 
URTO= 2. 302 Kurtosis 


9. At this point the calculator will automatically sort the 


data. This may take from several seconds to several minutes 


depending on the number of points in the data set. After 


the data set has been sorted, the calculator will display 
the median as follows: 


MED=2.400 TO (Press R/S) 
ee 2.500 


Two data values are displayed because when the number of 
data points is even, the median is not unigue, but rather 
spans an interval from the one point listed above to the 


other. Many users may wish to simply take the middle of 


this interval as the median, but any point is technically 
correct in the interval. When the number of data points is 
odd, the median is unique and only one value will be , 


displayed by the calculator. 


10. After the median is displayed as described in step 9, 
the calculator will display the following statistics labeled 


as shown; 


Display Meaning 
MINz2.100 Miniaum Value 
MAX22.900 Maxinugm value 
RNG20.800 Range 


11. At this point the calculator will respond with "CELL?" 


asking for the nusber of celis the user desires in the 
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histogram. If the number of cells is not significant at 
this point, the calculator will pick an appropriate number 


if the user simply presses run/stop. For this example, 


press: 
R/S 


12. Next the calculator responds with "WIDTH" asking for 
the width of the cells. Simply press run/stop if you do not 
wish to establish the width manually. Again, you may see 
the width the calculator will use by pressing the clear 
arrow key (Unless the width is an integer, you will also 
need to press FIX 3 to display the decimal properly if you 


wish to examine the width.) For this example, press: 


R/S 


13. The calculator will now display the cell frequency 
counts as an integer count followed by the next cell bound- 
ary. The leftmost cell boundary is set equal to the nininun 
value and is not explicitly output. If a data point falls 


exactly on a cell boundary, it is counted in the left cell. 


For this 2xample, the display 


will show: 


iii 


Display Meaning 
CNT=2 Two observations 
between 2.1 (the mi mm) 
XX=2.260 and 2.26 (the cell Sere een 
CNT=3 Three observations 
between 2.26 (see above) 
xx=2.420 and 2.42 (the next boundary) 
CNT=1 One observation 
between 2.42 (see abo s), 
xx=2.580 and 2.58 (the next boundary) 
CNT=3 Three observations 
between 2.58 (see a eh 
xx=2.740 and 2.74 (the next Bone ary) 
CNT=1 One observation 
between 2.74 (see above) 
xx=2.900 and 2.90 (the naxinum) 


14, After the last cell boundary is displayed, the calcuia- 
tor will display "CMD" asking the user for the next command. 
Prequently, the user will wish to modify the histogram by 

changing the number of cells or the cell width. To recalcu- 
late the histogram cell frequencies without re-entering the 


data press: 


XEQ ALPHA AGAIN ALPHA 


If no further work with the histogram is desired, the user 


may view the order statistics simply by pressing run/stop. 
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SINGLE VARIABLE STATISTICS 


DATA ENTRY LOOP 


(INCREMENT DATA STORAGE POINTER 
(RECALL DATA VALUE 


PREV. ENTERED? 


(PROMPT USER FOR NEXT DATA VALUE 


STORE THE DATA VALUE 
NO VERIFICATION OF DATA DESIRED? 


(FOLLOWING IS THE VERIFICATION ROUTINE 
(RECALL THE LABEL 


RECALL THE STORED DATA 

CLEAR DATA ENTRY FLAG 

WILL STOP POR DATA ENTRY IF F21 SET 
WAS THERE NO DATA CHANGE DURING VIEW? 
IF THERE WAS A NEW VALUE, THEN RECORD 
IT AND ACK AND RE-VERIFY THE DATA. 
POLLOWING IS THE STATISTICAL ACCUM. 
STORES SIGMA X 


(STORES SIGMA X-SQUARED 
(STORES SIGMA X-CUBED 
STORES SIGMA X-FOURTH-POWER 


IS THIS A DATA REVIEW? 
(IF DATA ENTRY, INCREMENT INPUT CNTR. 


(AT_END OF DATA ENTRY, RECALL INPUT 
COUNTER, WHICH IS A’ NUMBER EQUAL ONE 
MORE THAN 


NOMBR POINTS 
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HP41C SOURCE CODE: SINGLE VARIABLE STATISTICS 
( wee eee ow wee woe wwe ow oe ow 
"AGAIN DISPLAY HISTOGRAM 
, 
i I eres ne eee 
i 176 RCL 04 (ADDRESS BASE REGISTER 
173 
: + 
179 RCL IND x (RECALL THE FIRST ORDER STAT 
181 x29 97 CALL AN OUTPUT LABELING ROUTINE 
182 STO 09 HOLDS STARTING (LEFTHOS?) X BOUNDARY 
183 RCL 04 ADDRESS BASE REGISTER 
184 RCL 15 NUMBER OF DATA POINTS 
+ 
186 RCL IND x (RECALL THE N-TH ORDER STATISTIC 
188 XEQ 97 DISPLAY THE MAX VALUE OBSERVED 
189 RCL Z {NIN 
131 STO 08 (TEMP STORE POR THE RANGE 
193 XEQ 97 - (DISPLAY THE RANGE 
194 CP 00 INITIALIZE TEMP FLAG TO MARK LAST BAR 
135 RCL 15 NUMBER DATA POINTS 
196 RCL 04 ADDRESS BASE REGISTER 
198 1E3 
199 / jeonputite INDEX LOOP COUNTER 
290 RCL 04 ADDRESS BASE REGISTER 
| 303 ¢ 
i + 
204 STO 01 RO1 SET TO ADDRESS AND LOOP THRU DATA 
205 RCL 15 NOMBER POINTS 
25 
208 . ; (DEFAULT NUMBER OF BARS IS 2*LN(N) 
219 RND (VALUE IS ROUNDED NOT TRONCATED 
212 FROMPT {USER HAS OPTION TO CHNGE NOMBR CELLS 
213 RCL 08 RANGE 
214 X<>¥ 
ar “wrIDTH 
217 PROMPT USER HAS OPTION TO CHANGE CELL WIDTH 
218 STO 08 NOW HOLDS CELL WIDTH NOT RANG? 
333 Ber 08 CELL WIDTH 
221 ST + 09 OPPER LIMIT OF CURRENT CELL COUNTED 
234 STO 02 (INITIALIZE CELL COUNTER 
332 RCL IND 01 NEXT DATA POINT 
328 RCL 09 CELL UPPER LIMIT 
27 PIX 3 
228 RND 
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APPENDIX 3 


LINEAR PROGRAMMING EXAMPLE 


INTRODUCTION: 

Linear programming is an operations research technique 
normally associated with computerized data bases and the 
largest computers. Because of the complexity of the com- 
puter programs for linear programming and the large amount 
of data associated with most real world problems, calcula- 
tors have not been widely used for this application. With 
the increased memory capacity of the HP4ICV, however, it is 
now possible to offer a calculator program which can solve 
interesting small scale linear programs. Of value primarily 
as an educational aid, this program will also be able to 
solve many small scale problems found at Army division, bri-~ 
gade and battalion level. The text by Hillier and Lieberman 
{Ref. 8: pp. 16-66] is recommended as an introduction to 
the theory of linear programming as used by the program 
given in this appendix. Use of the program requires the 
user to formulate the linear programming problem; set up a 
Simplex tableau in standard form including adding slack, 


surplus and artificial variables as required; and interpret 
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the final tableau including the calculation of the values 
associated with the variables in the final basis. Using the 
tableau form of the Simplex algorithm, the calculator per- 
forms both phase I (to obtain a feasible solution) and phase 
II (to obtain an optimal solution) to solve the linear pro- 
gramming problem. The calculator automatically determines 
the pivot column and pivot row for each pivot step. Infeas- 
ible and unbounded problems are automatically idensified for 
the user by the programa. There is no explicit handling of 


variables with upper bounds. 


PROGRAM DESCRIPTION: 

The program is written as a series of subroutines, each 
of which performs a major step in the Simplex algorithm. To 
provide clarity to the user, alphabetic labels have been 
retained to identify the subroutines in lieu of faster and 
more memory efficient numeric labels. The alphabetic labels 
have not been retained for use as program entry points and 
may be changed to numeric labels at the option of the user. 
The program has two entry points, "LP" for running a new 
problems and "ALP" for reviewing data previously entered. 

Subroutine "FPINDQ" determines the pivot coluagn by 
selecting the variable in the objactive function with the 


gost negative "price." If "PINDQ" discovers at least one 


negative value in the objective function, then the tableau 
column number associated with the most negative value will 
be stored in register 05. Upon return from "PINDQ,” the 
main routine tests register 05 to see if it contains a non- 
zero entry. If the entry is zero, it means that no further 
pivots will improve the value o% the objective function, and 
the Simplex algorithm halts. If the program was in phase I 
(flag 11 clear) and the value of the phase I objective fuac- 
tion is reduced to zero, then a feasible solution has been 
found and the program will automatically proceed to phase II 
to discover an optimal solution. 

Subroutine “"PINDP"” determines which variable will leave 
the current basis by performing a minimum positive ratio 
test along the pivot column. In this way, the pivot row is 
determined. The row number of the pivot row is stored in 
register 06. Upon return from subroutine "FINDP," the main 
routine tests register 06 to see if it contains a non-zero 
entry. If the entry is zero, it aeans that the problerp is 
unkounded and the Simplex algorithm halts. Such an 
unbounded condition is most likely caused by an error in the 
probles formulation. 

Having determined the pivot column and the pivot row, 


subroutine "PIVOT" performs the actual Simplex pivot 
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operation. To speed calculation register 00 is used as a 
temporary register to hold the reciprocal of the pivot 
element. Note that the pivot row is handled separately from 
the other rows in the tableau. lag 04 .s used to provide 
the option of stopping calculation after every pivot. When 
this flag is set, the program will halt to allow the user to 
review the status of the tableau with the "ALP" function. 
Subroutine "CHECK" has two primary functions. First, it 
is used to verify that the designated basic variables are in 


crow elimination form prior to the start of the Simplex 


i ae ce ee ane a 


algorithm. This means that the basic variable must have a 
coefficient of 1 in the row in which it is basic and zero's 
in all other rows. The second function of check is to pre- 
pare the objective function for phase I, if the initial 
basis contains artificial variables as indicated by one or 
more minus signs in the "JB" vector. 

Three subroutines are used to query the user for input 
data. Subroutine "READMN" queries the user for the number 
of constraints and decision variables in the problem and 
verifies the calculator memory is set to contain all the 
data necessary to solve the problea. Subroutine "READJB" 
queries the user for a column vector of pointers which indi- 


cate which variable is currently basic in each row. When 
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entering this vector of pointers, the user indicates 
artificial variables with a2 minus sign. Subroutine “READA" 
queries the user for the values in the initial Simplex 
tableau including the slack and surplus variables and the 
right hand side and objective function. 

Several other service routines also are provided in this 
program. Memory size verification is done by subroutine 
"SIZE," which is called from within "READMN." Subroutine 
"IN" is used to query the user for data entry and is called 
by all of the data input routines. Subroutine "NXT" ini- 
tializes registers which contain frequently used quantities 
such as the the total size of the tableau for phase I and 
phase II. Subroutine "INIT" clears the calculator memory 
and sets flags and program constants appropriately for input 
of a new problem. Subroutine "SETL" establishes the loop 
counters used repeatedly within almost every other subrou- 
tine. Subroutine "ERR1" displays an appropriate error 


message if a data entry error is detected. 


SANSPLE PROBLEM: 


A division assistant G4 is planning an ammunition upload 


plan. There are four types of tank aunitions to consider, i 4 


including: 
A = Discarding Sabot Rounds 
B = High Explosive Anti-Tank Rounds 
C = Phosphorous Munitions 
D = Machine Gun Ammunition 


Based on the Commander's guidance the assistant G4 is to 
consider the sabot rounds as twice as important tne HEAT 
rounds, which in turn are themselves twice as important as a 
unit amount of phosphorous sunitions and machine gun 


amaunition. His mission then, is to maximize: 
Z= WA +e 2B + C+D 
He is, hewever, constrained by the following factors: 


1. There can be no more than 30 units of both sahbot 
and HEAT munitions combined. 


2. There can be no more than 50 units of all types of 
amaunition combined. 


3. There must be at least 30 units of HEAT and 
phosphorous munitions combined. 


4. There must be at least 5 units of machine gun 
agasunition. 


These constraints may be expressed as; 


A+ B s 30 
A+ BeCe#DdD € 50 
Bec 2 30 

D 2 5 


Based om the Commander's guidance and the constraints listed 


above, formulate an optimum load plan. Practional units are 


permitted. 


SOLUTION: 

1. Before beginning with the calculator, the first step is 
to layout the tableau in standard form. This step and the 
last step of interpreting the final tableau require working 
knowledge of linear programming as explained in Hillier and 
Liberman {[Ref. 8: pp. 16-66jJ. The standard form of the 


tableau is: 


1 2 3 4 5 6 7 8 9 10 11 
A 3 c D H1 H2 $1 S2 Al A2 RHS 

1 1 30 

1 1 1 1 1 50 

1 1 -1 1 30 

1 -1 1 5 


Ses Sa Dy Sea aE aaa ad a a ead ¢ (aaa! 
In this tableau, H1 and H2 are surplus variables; $1 and S2 


are slack variables; and Al and A2 are artificial variables. 


2. The first step with the calculator is to set the size of 
the calculator's data memory. This program requires 20 reg- 
isters for teaporary storage, 1 register for each tableau 
element, and 1 register for each row to hold the pointer to 
the basic variable for that tow. Thus, if MN is the number 
of constraints and N is the number of variables including 


slack, surplus and artificial variables, then the total data 


Te ee ee 


storage requirement is: 
storage required = 21 + M + ((N + 1) x (M # 2)) 


As mentioned in the program description, the "SIZE" subrou- 

tine will automatically verify that the user has allocated 

enough data storage to solve the problem. The length of the i 
program is such that 177 data storage registers is the maxi- 
mum number of data storage registers that can be allocated. 
Thus, linear programs with 7 constraints and 15 decision 
variables can be solved with this program. For this 


example, press: 


XZQ ALPHA SIZE ALPHA 175 


3. Call for execution of the program with a new data set. 


Press; 


XEQ ALPHA LP ALPHA 


4, The calculator will respond with the prompt "NUM ROWS?" 
asking for the number of constraints in the linear program 


foraulation. In this example, press: 


4 R/S 


5. The calculator will respond with the prompt "NUM COLS?" 


asking for the number of variables in the problem. The user 


must count the number of slack, surplus and artificial 


variables in this total. In this example, press: 
10 R/S 


6. The calculator will respond with the prompt "BASIC 1 2" 
asking for the variable number of the variable which is 
basic in the first row. One of the major features of this 
program is that the basic variables need not be in the 
rightmost positions in the tableau. Thus, if a tabieau were 
given in which some pivots had already been performed, the 
program could resume operation immediately. In this 


example, press: ~ 
7 RSS 


In a Similar fashion, the calculator will then query the 
user for the variable number of the variables which are 
basic in the remaining rows. 


For this example: 


see Respond with 
Basic 3 2 8 R/S 
Basic 2? 9 CHS R/S 
Basic 4 ? 10 CHS R/S 


Notice that because the basic variables in rows three and 
four are artificial variables, the variable number is 


entered as a negative number. This signals the calculator 


that these variables must be driven from the basis in order 


to reach an initial feasible solution. 


7. \Next, the calculator will respond with "T1,1?" asking 
for the first element in the tableau. The usar should enter 
the numbers in the tableau using the digit entry keys and 
pressing run/stop after every entry. Notice that the right 
hand side and the objective function will be entered with 
the appropriate index in the tableau as shown in step 1 
above. The user must insure that the objective function is 
in standard form with the appropriate sign for each 


coefficient--in this example each coefficient is negative. 


8. After the last element in the tableau has been entered, 
the calculator will begin to automatically perform the Sin- 
plex algorithm. If the user wishes to stop the calculator 


after every pivot, he may at any time press: 


R/S SP 04 R/S 


If this flag is set, the calculator will stop and display 


the pivot number after every pivot is completed. 


9. When the Simplex algorithm can no longer improve «he 
objective function, the calculator will stop and display the 


value of the objective function. In this example, the 
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cng tem tena tens 


calculator will stop after approximately three minutes and 


display: 
VALUE=110.000 


10. At this point, the user must use entry point "ALP" to 
determine the status of the final tableau. For this 


example, press: 
XEQ ALPHA ALP ALPHA 


Then by sequentially pressing the run/stop and clear arrow 
keys, the basic variables and tableau entries will be 


displayed. Por example, in this problen: 


See Press see Meaning 

BASIC 1? CLX 2 Variable 2 is basic 
in the first row. 

BASIC 2? CLX 1 Variable 1 is basic 


in the second row. 
etc. 


Then for the elements of the tableau: 


see Press see Meaning 
T1,1? CLE 0.000 Tableau entry 
T1,2? CLX 1.000 Tableau entry 


11. After the calculator is finished, it remains for the 
user to interpret the final tableau. Again, the reference 
by Hillier and Lieberman [Ref. 8: pp. 16-66] is of primary 
value. In particular, the user must be able +o determine 


the value of the final decision variables based upon what 
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variables are in the basis, and what the final "right hand 


seen ari 2 ae On 


tableau is: 


-1 


1 


side" values are for each row. 


“1 


For this example, the final 


7 48 9 10 

$1 $2 Al A2 

1 -1 1 
1-1-1 
1 -1 


1 


11 
RHS 


15 
15 
15 

5 
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Thus, the solution may be interpreted as 15 units each for 


munitions A,B and C and 5 units for aunition D. 
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INTRODUCTION: 


The calculator subroutines described in this appendix 
perform functions which are frequently required by applica- 
tion programs and are therefore ideal candidates for use in 
a read only memory (ROM.) These routines were written espe- 
cially to illustrate the differences between read only mea- 
ory routines and routines designed for individual use via 
bar code or magnetic cards. These differences include more 
attention to entry and exit point options, an attempt to 
keep the size of the routines as compact as possible, and an 
attempt not to disturb the register stack if at all 
possible. 

These common subroutines are provided separately from 
application programs because when more than one application 
program uses the routines, as is recommended for these pro- 
grams, the use of a separate block of common functions saves 
space in the ROM overall. Also, by providing a convenient 
set of "sgacro" instructions, application pregrams can be 


constructed more guickly and easily. Because these 
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subroutinas are used frequently, they have been individualiy 
optimized and tested to save memory space and executicn 
time. By using «hese “macros within an application pro- 
gram, the application programmer can be reasonably certain 
of their efficiency and reliability. 

Almost all user/calculator interface is handled by these 
routines. There is one set of subroutines which assumes the 
user has a printer, and one set which does not. Printer 
instructions are preceeded in the listings shown in the 
appendix by an (PRT: label. When not using these routines 
on read only memory, the user will load one set or «he other 
(but not both), as appropriate for his/her application. In 
so doing, the user with the printer gats full benefit from 
it while the user without the printer pays no penalty in 
execution time or memory space for the calculator's print 
instructions. Also, to change from use of the printer to 
use of the calculator without it, the user needs only to 
read in the new common block~-the application programs are 
retained in memory unchanged. The subroutines appear «he 
same to any application program--giving the added benefi+ 
that any application program which uses them for input or 


output operations will automatically make good use of 


the printer even if written by a programmer who did net 
\ explicitly consider a printer requirement. 

When using these common subroutines, a discipline is 
enforced upcn the application program concerning use of the 
calculator memory registers. This saves the programmer 
from having to plan his “register map" from scratch for each 
new program. Also, it insures compatability across differ- 
ent application programs for similar data objects such as 
matrices and loop index counters. One of the most annoying 
probleas with read only memory programs available from the 
calculator manufacturer is this lack of cross program com- 
patability. Conflict in the use of memory registers is the 
rule, rather than the exception; and it is frequently impos- 
sible to efficiently use more than one read only memory pro- 
gram as a subroutine within a user written program. A ‘third 
reason why register assignment standards are advantageous is 
that they make it easier for the user of the calculator to 
remember the key register assignments and, if necessary, 

recall their contents during the execution cf an application 
prograa. 
Another function performed by this set cf common subrou- 
tines is to simplify the use of indirect addr2ssing--a 


critical goal on the AP4ICY. 


Because the comaon subroutines Listed in this appendix 
are always called by application programs and never from the 
keyboard by the user, the typical user instructions are 
inappropriate. Instead, for the benefixc of application 
programgers wishing to use the routines, the basic functions 
and structure of each are explained in subsequent sections 


of this appendix. 


subroutine "IN" 

Subroutine “IN® is used as a general input and output 
interface between the user and the calculator. This subrou- 
tine has three alternative entry points which when called 
affect functions as follows: 


eas mode (displays a ignest tion mark guery) 
O--Output mode (displays labeled data value) 
IX--Direct node inpat Of value in x register) 


In particular, one entry point, "IN", nay be called whenever 
an application program must guery the user for a numeric 
input value. as such, it is a direct replacement for the 
PROMPT instruction organic to the calculator, but automati- 
cally proupts, verifies and stores the received value using 
an indirect address contained in register 00. The printer 
version of the subroutine will automatically label and print 


the final, verified data value recorded. 
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Subroutine "IN" uses register 00 as a data location 
pointer and automatically increments this register so that 
subsequent calls to the subroutine will resul+ in sequential 
data manipulation. The application programmer must insure 
that register 00 contains a number equal to the storage reg- 
ister number prior to calling the subroutine. For example, 
if register 00 contains 17, “IN" will store the data in reg- 
ister 17. One of the major advantages of this esuv ine is 
that the same subroutine may be used to verify and/or change 
the data previously recorded. Thus, separate edit routines 
are usually unnecessary. Pressing the R/S key without 
touching any other key leaves the value stored unchanged. 
Pressing "1" and *¢* and then "R/S" adds one to the current 
stored value, and so on for other functicn keys. Entering a 
new string of digits results in that new value being stored. 

An additional feature of this subroutine is the "verify" 
mode indicated by flag 05. Flag 05 is reserved for this 
purpose and is set "on" by a call to subroutine "VRN-- 
another of the subroutines in this common set. The verify 
mode is intended for use by a novice or Sekar user who 
wishes to verify every data value as it is keyed into the 
calculator. The advantage is increased accuracy and 


confidenca in the result. 
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Subroutine "D2" 

Subroutine "D2" is used to set up the indax register for 
a program loop. This subroutine has two alternative entry 
points which when called, increment different index regis- 
F ters as follows: 


tablishes register 02 as the index 
tablishes register 67 as tha index 


This subroutine is intended for use with the "ISG" loop 
structure which has the effect most like that of a FORTRAN 


"DO LOOP." Por exaaple, to execute a loop 20 times: 


The advantage of this form of loop structure is that 


eS 


register 00 may be used within the loop for address calcuia- 
tions. The first time the loop is addressed the integer 


portion of the number in register 00 will be 1, the second 


eet 


time it will be 2 and so on. There is no need to «runcate2 


the fractional portion of the number because the HP41C 


ignores the fractional component of a number when 
calculating a register address. JUse of index registers for _ P 


address calculation makes indirect addressing practical. 


Sat aenaldee 
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Registers 01 and 02 should be reserved for use as indey 
tegisters by the application programmer. In most cases 


these two registers should prove sufficient. 


Subroutine "YR" 

Subroutine "VR" is used as a general purpose calculator 
initialization routine. This subroutine has «hree al+erna- 
tive entry points which vary the amount of initialization 
performed as follows; 


VR~-Sets the verify mode on, and the following; 
wR--Suppresses the audio tones, and the following: 
WS--Clears all measory : 

Sets the displa for integers, 

Assigns statistical regisers . 

Sets "zero" level for equality testing, and 

Sets base address for indirec addrassing. 


In the printer version of this initializatior routine, 
the subroutine prints a banner (usually the program name) 
which has been stored in the alpha register prior to calling 
the initialization subroutine. 

If flag 13 is set prior +o calling the initialization 
routine, <hen the calling program must have placed the nua- 
ber of data registers required to execute the program in the 
x register prior to calling the initialization routine. In 
this case, a check will automatically be performed using 
subroutine "SZ" described below. 

It is recommended that all application programs provide 


an alternate entry point which bypasses the initialization 


step. Then if a data error is encountered, the user may 
review the data entered into the calculator by simply press- 
ing the return key Once after every prompt. This procedure 
works because subroutine "IN" recalls the stored value prior 
to proapting the user. When the user presses the clear key, 
the alphabetic proipt is removed and the existing data value 


revealed. 


Subroutine "$2" 

Subroutine "SZ" is used to test if sufficient numbers of 
data registers are available ¢o run an application prograg. 
This subroutine may be either called directly or as part of 


the initialization routine described above. 


Subroutine "ER" 

Subroutine "ER" is called whenever the application pro- 
gram encounters an error--usually in the input data. A 
prompt is displayed and an audio tone sounded, provided flag 
26 (the audio enable flag) has not been cleared by the 


initialization routine described in paragraph D. 


aubroutine "SORT" 


This subroutine is included to illustrate the use of a 
stack register table in the program comments, but it also 


represents a useful utility routine. The sorting algoritha 


a 
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used is the shell sort [Ref. 6: pp. 84-95] which gives 
reasonably fast sorting times with a very small program 
size. All conventions such as base register in ROY and num- 
ber of data points in R15 are assumed by this subroutine. A 
cosplete list of all such register assignments is listed a+ 


the end of the program listing. 


Subroutines "PUI" and "GET" 

These two small routines provide a useful capability to 
store and recall up to three integers between 0 and 999 in 
one data register. This means that if you ars manipulating 
a spread sheet of small, positive integers you can store the 
same data in one third the space. Of course, run time is 
degraded (about 20 seconds for every 100 data references.) 
To store a value, assuming the base register has been 


defined, just press: 
value ENTER| point-number XEQ "PUT 
To recall a value, siuply press: 


point-nuaber XEQ "GET 
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APPENDIX D 
THE CROSS COMPILER PROGRAN AND COMMAND PROCESSOR 


DESIGN METHODOLOGY: 

This appendix discusses the design methcdology used dur- 
ing construction of both the cross compiler program and the 
command processor, which is an IBM EXEC II program which 
provides an interactive programming environment for users of 
the systen. 

Blazie's compiler for the HP65 calculator {Ref. 9] rep- 
resents one of the first attempts to provide a compiler for 
calculator programs. Both Carvalho [Ref. 10: pp. 25-29] 
and McNeal (Ref. 117: pp. 148-178) have published BASIC lan- 

| guage programs which cross compile HP4ICY instructions on a 

| microcomputer for output to a line printer which can print 
acceptable kar code. While these referenced programs pro- 

) vided valuable insights into the problem, especially into 
the special characteristics of the HP4ICV instruction set, 
none was 2xactly suited to the needs of this study. Pecause 
the Versatec plotter at the Naval Postgraduate School could 
be easily used only by FORTRAN programs, FORTRAN seemed the 


coaputer language of choice for this project. Both prograas 
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were written with limited objectives and neither would have 
aasily suppcerted the extensions desired. Extensions planned 


for implementation included: 

e An extended instruction set. 

e In code comments. 

e Extensive error checking. 

e Compatability with the 2aulator. 
® Synthetic Instructions [ Ref. 1}. 


e Instruction macro's. 

Having decided to code an original cross compiler, a 
design methodology which would capitalize on «he advantages 
of FORTRAN was planned. FORTRAN'’s major deficiancy for use 

| in constructing a compiler of any type is its lack of 
alpha-numeric string handling capabilities. Rather than 
struggle with the lack of string functions, it was decided 
to code the necessary string functions as separate subrovu- 
tines. This decision was reinforced countless tines 
throughout the process of writing the compiler. The string 
function subroutines have been used not only in the cross 


compiler, but in many other FORTRAN prograas since they were 


originally written. In fact, many persons who have no 


interest in the HPUICY cross compiler may find the set of 
string functions listed in this thesis to be a valuable set 
of utility routines to be used +o augment FORTRAN. The gen- 
eval convention used througout the string function subrou- 
tines is that an alphabetic string may be represented as a 
vector of two byte integer variables used to store the char-~ 


acters and a single four byte integer variable used +o 


2 Re CE eee ee 


store the length of the string. 
One of the major advantages of the cross compiler is its 
ability to handle comments integrated within the HP4ICV 
source code. This feature is critical to the clarification 
of the logical structure of the HP41CV programs. Because 
the parenthesis is not a valid HP41CV character, it was cho- 
sen as tha comment indicator character. A comment may occur 
beginning at the first column on an input line or anywhere 
after an HP41CVY instruction. The comment must follow the 
instruction hecanse everything after the comment mark out to 
the eng of the input line is considered part of the ccmment. 


The control the user has over the output listing is also 


; 
' 
{ 


ment indicator marks are placed in positions one and two of 
i 


one of the advantages of the cross compiler. When ¢we con- 


the input line, the compiler will force a page eject when 


printing the output listing. In addition, the user can vary 


the number of output lines per page and cause useful kLanners 
to be placed adjacent to program labels for ease of 
recognition. 

Altogether there are twenty-four subroutines and a main 
program which constitute the cross compiler. The source 


code for each of these routines is provided in the second 


a ee 


section of this appendix. Each subroutine begins with a 
statement concerning its function and construction. Accord- 
ingly, no general description of each subroutine will be 
repeated here. However, subroutine COMP deserves special 
attention, for it is the master lexicographic analyzer for 
the compiler and would also interface the user with the eau- 
lator. Its function is to receive a single line of HP4 ICV 
source code and identify it. COMP considers all HP4&ICV 
instructions to be of one of three types. The firs*+ cat- 
egory are the single byte instructions with no operands that 
can be compiled by a simple table look up. COMP has been 
constructed so that the instruction set can be extended at 
any time simply by increasing the size of this table. [hn 
this way abbreviated or altered command names could be eas- 
ily used. The second category of instructions are the 


multi-byte instructions which require a table lookup and the 


translation of one or more operands, including possibly an 


indirect instruction indicator. fhe table examined by thea 
compiler is the same as for the category one instructions, 
and a code is given in the table which indicates to the 
compiler the number of operands which are required with each 
instruction. A syntax check is then made in subroutines 


IONE and ITWO to insure that the number and characteris- 


tics of each operand are appropriate for the given instruc- 
tion. One of the major advantages of the use of the cross 
compiler is the syntax and error checking that is performed 
during the compilation process. fhe third type of instruc- 
tion reprasents the exceptional instructions that are so 
difficult to compile that they reguire separate subroutines 
for efficient coapilation. These instructions include stor- 
age and recall of data, program labels and program flew con- 
trol statements such as goto and axecute. 

In order to provide an efficient prograaming command 
system for the compiler that would minimize the need to know 
technical details abou the operation of the compiler, an 
Iam EXEC II program was written. This program not only 
interfaces the user to the compiler, but it also provides on 
line user instructions as to how to use the systen. 


Included in this command processor is a command menu which 


gives the format and short description for each command. 


Tose AMA et ne 6 mpg ol 


Another command, help, provides more detailed information 
about each command. When a novice user first enters the 

exec, or types the name of the exec followed by a question 
mark, then he receives a four page narative description of 


what the system is, how it works, and what actions he must 


take to write a successful HPGICV application progran. 
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